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Abstract 

The contribution of a spin-rich qqq force ( in conjunction with 
pairwise qq forces) to the analytical structure of the qqq wave func- 
tion is worked out in the high momentum regime of QCD where the 
confining interaction may be ignored, so that the dominant effect is 
Coulombic. A distinctive feature of this study is that the spin-rich 
qqq force is generated by a ggg vertex ( a genuine part of the QCD 
Lagrangian ) wherein the 3 radiating gluon lines end on as many quark 
lines, giving rise to a (Mercedes-Benz type) Y-shaped diagram. The 
dynamics is that of a Salpeter-like equation ( 3D support for the ker- 
nel) formulated covariantly on the light front, a la Markov- Yukawa 
Transversality Principle (MYTP) which warrants a 2-way intercon- 
nection between the 3D and 4D Bethe-Salpeter (BSE) forms for 2 
as well as 3 fermion quarks. With these ingredients, the differential 
equation for the 3D wave function <f> receives well-defined contribu- 
tions from the qq and qqq forces. In particular a negative eigenvalue 
of the spin operator iu\.a-i x 03 which is an integral part of the qqq 
force, causes a characteristic singularity in the differential equation, 
signalling the dynamical effect of a spin-rich qqq force not yet consid- 
ered in the literature. The potentially crucial role of this interesting 
effect vis-a-vis the so-called 'spin anomaly' of the proton, is a subject 
of considerable physical interest. 
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1 Introduction 



The concept of a fundamental 3-body force (on par with a 2-body force) 
is hard to realize in physics, leaving aside certain ad hoc representations of 
higher order effects, for example those of A, N* resonances in hadron physics. 
At the deeper quark-gluon level on the other hand, a truly 3-body qqq force 
shows up as a folding of a ggg vertex ( a genuine part of the gluon Lagrangian 
in QCD) with 3 distinct qgq vertices, so as to form a F-shaped diagram (see 
fig 1 below). Indeed a 3-body qqq force of this type, albeit for 'scalar' gluons, 
was first suggested by Ernest Ma when QCD was still in its infancy. [ A 
similar representation is also possible for NNN interaction via ppp or aaa 
vertices, but was never in fashion in the literature [2]]. We note in passing 
that a F-shaped (Mercedes-Benz type) picture [3] was once considered in the 
context of a preon model for quarks and leptons. 

In the context of QCD as a Yang-Mills field, a ggg vertex has a momentum 
representation of the form [I] 

Wggg = -Wsfabc[{h ~ k 2 ) x 5^ + (h - k 3 )^5 uX + {h ~ h)Jxfi] (1-1) 

where the 4-momenta emanating from the ggg vertex satisfy k% + k 2 + k$ = 0, 
and f a b c is the color factor. When this vertex is folded into 3 qgq vertices 
of the respective forms g s u(p' 1 )ij^{X1/2}u(pi), and two similar terms, the 
resultant qqq interaction matrix (suppressing the Dirac spinors for the 3 
quarks) becomes 

V m = f|[*7 (1) .(^ - ^ 3 )7 (2) .7 (3) + {2} + miMXM/ikfklkl} (1.2) 

where k{ = Pi — p'^ A« are the color matrices which get contracted into the 
corresponding scalar triple products in an obvious notation. [Note that the 
flavour indices are absent here since the quark gluon interaction is flavour 
blind]. 

This interaction will be considered in conjunction with 3 pairs of qq forces 
within the framework of a Bethe-Salpeter type dynamics to be specified be- 
low. Before proceeding in this direction, it is in order to explain a possi- 
ble motivation behind the use of a direct qqq force with such a rich spin 
dependence. Apart from the intrinsic beauty of this term, the immediate 
provocation for its use comes from the issue of "proton spin" which, after 
making headlines about two decades ago, has come to the fore once again, 
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thanks to the progress of experimental techniques in polarized deep inelastic 
scattering off polarized protons, and their variations thereof, which allow for 
an experimental determination of certain key QCD parameters by relating 
them to certain observable quantities emanating from external probes ; ( see 
a recent review [5] for references and other details). On the other hand it 
is also of considerable theoretical interest to determine these very quantities 
directly from the intrinsic premises of QCD provided one has a " good " 
qqq wave function to play with. Such a plea would have sounded rather 
Utopian in the early days of QCD when phenomenology was the order of 
the day. Today however many aspects of QCD are understood well enough 
to make such studies worthwhile by hindsight, with possible ramifications 
beyond their educational value. The role of the direct 3-body force may be 
seen in this light, while working for simplicity in the experimentally accessi- 
ble regime of valence quarks. To that end this paper is specifically concerned 
with the effect of the 3-body force (1.2) on the analytical structure of the 
qqq wave function, while the formalism dealing with the contributions of the 
various operators ( 27^75, 2-gluon effects, etc) to the proton spin is reserved 
for a subsequent paper. 

1.1 Theoretical Ingredients 

In the valence quark regime, we need to consider a qqq system governed by 
pairwise qq forces as well as a direct 3-quark force of the type (1.2). A further 
simplification occurs in the high momentum regime where the effect of confin- 
ing forces may be neglected, so that only coulombic forces need be kept track 
of. We shall take the dynamics of a qqq hadron in the high — momentum 
regime to be governed by a Bethe-Salpeter Equation (BSE) whose kernel 
is a sum of three pairs of (coulombic) qq forces plus a single 3-body term 
(1.2). Unfortunately the 4D form of the BSE is too general to be of prac- 
tical value for qqq dynamics, so a better option is the Salpeter Equation [6J 
representing its instantaneous form. And, except for its lack of covariance, 
the Salpeter equation has the remarkable property of 3D-4D interlinkage, a 
feature that had been present all along in the original formulation itself [BJ, 
but had somehow remained hidden from view in the literature, until clarified 
[7] in the context of a comprehensive two — tier BSE formalism developed 
independently in a covariant manner [Sj. (A practical significance of the 
'two-tier' formalism is that the reduced 3D form is ideal for the determi- 
nation of hadron mass spectra of both qq and qqq types [HI [TD|, while the 
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reconstructed 4D form is convenient for the evaluation of transition ampli- 
tudes [5] via Feynman techniques for loop diagrams; see also Munz et al 
ref.(9) ). A covariant formulation of the BSE is centered around the hadron 
4-momentum P M in accordance with the Markov- Yukawa Transversality Prin- 
ciple (MYTP) [Til [12] . which was shown to be a 'gauge principle' in disguise 
|13j . It ensures that the interactions among the constituents be transverse 
to the direction of P^. In the high momentum regime to be considered here, 
the confining interaction has been ignored for simplicity, which leaves the 
3D form of the BS dynamics inadequate for mass spectral determination, yet 
its dynamical on the spin-structure of the wave function should be realistic 
enough for dealing with the hadron spin in the high momentum limit. 

Now to another vital element of the theory: Although the Salpeter Equa- 
tion, as the instant form of BSE, admits of a covariant formulation in the 
rest frame of the total hadron 4-momentum P^, it suffers from certain ill- 
defined 4D loop integrals due to a 'Lorentz-mismatch' among the rest-frames 
of the participating hadronic composites, resulting in time-like momentum 
components in the exponential/gaussian factors associated with their vertex 
functions. This is especially true of triangle loops, such as applicable to the 
pion form factor and p — tttt coupling [H] where this disease causes unwar- 
ranted "complexities" in the amplitudes, while two-quark loops just escape 
this pathology. For a possible remedy against this disease, without losing the 
benefit of an MYTP - based 3D-4D interconnection, a promising candidate 
is the light-front approach of Dirac [12] by virtue of its bigger (7) stabil- 
ity group compared with 6 for the instant form theory. Its basic simplicity 
was first noticed by Weinberg [16] who formulated the infinite momentum 
frame towards the same end. A covariant generalization of MYTP on the 
light front, requiring the use of two null 4- vectors n^ and that satisfy 
n 2 = n 2 = and n.n = 1, cures the 'Lorentz mismatch' disease noted above 
[j~7] . The remaining problem of n- dependence of transition amplitudes is 
solved through a simple ansatz of 'Lorentz completion' [HI [18], so as to yield 
a Lorentz-invariant pion e.m. form factor in accord with experiment . A sim- 
ilar approach had been considered by Carbonell et al [32] in the context of 
the Kadychevsky-Karmanov formalism [2U1 EE] , except for their missing out 
on the second (dual) null- vector h which happens to be crucial for recovering 
a Lorentz - invariant structure in a more natural way. 
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1.2 Plan of the Paper 

The plan of the paper is based on an interlinked 3D-4D BSE formalism char- 
acterized by a Lorentz-covariant 3D support for its kernel a la MYTP [TT| [T2] , 
adapted to the light front (LF) [T7] . The full structure of the kernel is a sum 
of 3 pairs of coulombic qq forces plus a qqq force, Eq.(1.2), whose 3D support 
is implied by the fact that all internal momenta q be transverse to P^ , 
viz., % = q^ — q.PP^/P 2 . However the propagators are left untouched in 
their standard 4D forms. [The light front formulation requires a collinear 
frame [17] which is further specified in Section 2 ]. The strategy now lies 
in a step-wise reduction of the (Salpeter-like) 4D Master Equation involving 
the actual (4D) fermionic wave function Step A consists in expressing 
\I/ in terms of an auxiliary (bosonic) 4D quantity $ satisfying an equivalent 
(bosonic) 4D (Salpeter-like) equation. Step B involves a " Gordon reduction 
" to eliminate the 7 M matrices in favour of a^ v matrices. Finally Step C con- 
sists in a 3D reduction of this 'bosonized' Master Equation, and a subsequent 
reconstruction of the original \I/ in 4D form by a suitable reversal of steps. 
In the process a 3D scalar wave function is introduced which not only fa- 
cilitates an explicit solution of the 3D (albeit fully covariant) Salpeter-like 
equation but is also a key component of the (reconstructed) 4D fermionic 
wave function To facilitate the process of 3D-4D interlinkage, a Green's 
function approach is employed a la [22], from which it is straightforward to 
derive the corresponding wave functions via the appropriate 'pole' limits. 
The entire exercise involves a close correspondence between the (earlier) in- 
stantaneous form [22] and the (later) LF form [T7J, so as to project only the 
latter by making free use of the results of the former. A new element will be 
a generalization of the earlier formalism [22] so as to include the effect of the 
3-body force, Eq (1.2), on the structure of the qqq Green's function which 
will require a more elaborate strategy, keeping in view the relative strengths 
of qq and qqq forces. 

After a short account of the correspondence between the instant and LF 
forms of the dynamics, Section 2 is devoted to Steps A and B for converting 
the Master Equation for the actual 4D qqq wave function \l/ (fermionic) to 
an equivalent 4D (bosonic) form involving $, with all 7- matrices eliminated 
in favour of Pauli matrices, by exploiting the effective 3D support for the 
qq and qqq kernels on the light-front, as described above. Section 3 de- 
scribes the first part of Step C, viz., the use of Green's- function method 
for the 3D reduction of the Master Equation for the 4D Green's function 
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G corresponding to $, resulting in an integral equation for the 3D Green's 
function G corresponding to phi, specialized to the large momentum regime 
on the light-front. The second part of Step C, namely its reconstruction 
back to the 4D quantity G a la (22], leading to a formal connection between 
$ and 0, is the subject of Section 4, now with the added complexity behind 
the inclusion of the 3-body qqq force along with the (dominant) pairwise qq 
forces. [ The subsequent reconstruction of the 4D fermionic form follows 
from the results of Section 2]. As a further aid to the understanding of the 
3D-4D interconnection, Appendix A describes the complete procedure for a 
prototype qq subsystem, whose results are freely used for justifying several 
steps in Sections 3 and 4. In Section 5, the full 3D BSE for <fi is set up in 
a simplified form designed to check for the dynamical effect of spin present 
in both the qq and qqq forces, on the structure of the differential equation 
for <p (the dominant effect being of the latter ! ) . Section 6 is devoted to 
a short critique of the role of the V qqq term on the analytic srtucture of the 
qqq wave function, while the derivation of the requisite QCD parameters for 
the baryon from the forward scattering amplitude off the ^7^75 operator, 
via the quark constituents, as well as a more elaborate 2-gluon contribution 
to the proton spin, is reserved for a subsequent paper. 

2 Salpeter Eq on LF : From ^ To $ 

2.1 Master Eq : Instant Vs LF Forms of Dynamics 

Since the first step in the formulation of a covariant Salpeter-like equation 
on the light front (LF) is to establish a correspondence between the instant 
and LF forms of the dynamics, we first recall some definitions [17] for the LF 
quantities p± = po=tP3 defined covariantly as p + = n.p^/2 and p_ = —n.p\/2. 
while the perpendicular components continue to be denoted by p± in both 
notations. Now for a typical internal momentum g M , the parallel component 
P.qP^/P 2 of the instant form translates in the LF form as = zPnU^, 
where P n = P.h, and z = n.q/n.P. As a check, q 2 = q\ + z 2 M 2 which shows 
that zM plays the role of the third component of q on LF. Next, we collect 
some of the more important definitions / results of the LF formalism [UJ 

q± = q ~ q n n;q = q± + zP n n; z = q.n/P.n;q n = q.h; (2.1) 
P n = P.n;P.q = P n q.n + P.nq n ;q.n = P±.q± = 0; 
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P.q 



P n q.n; q 2 = q\ + AfV; P 2 = -M 2 



For a qqq baryon, there are two internal momenta, each separately satisfying 
the relations (2.1). Note that for any 4- vector A, A.n and — A n correspond 
to 1/ \[2 times the usual light front quantities A± = A ± A z respectively. 
But since a physical amplitude must not depend on the orientation n, a 
simple device termed Lorentz — completion via the collinear trick [17] yields 
a Lorentz-invariant amplitude for a transition process with three external 
lines P, P', P"{= P + P') as follows. Since collinearity implies P±.P[_ = 0, the 
4-scalar product P.P' = P.nP^+P n P' .n+P ± .P± simplifies to P.nP^+P n P'.n. 
Then ' Lorentz completion ' simply amounts to reversing the last step via the 
'zero' quantity P±.P' ± , so as to recover the Lorentz invariant quantity P.P' at 
the end ! And as a practical simplification, one does not even have to use the 
n, n symbols; it suffices to use the more familiar light-front components A± 
for the covariant quantities v2[A.n, —A n ] respectively. For ready reference, 
the precise correspondence between the instant and LF definitions of the 
'parallel (z)' and 'time-like (0) ' components of the various 4- momenta for 
a qqq baryon ( i = 1,2,3) [221 ES] : 



The last part of Eq.(2.2) defines a covariant 3- vector on the LF that will 
frequently appear in the reduced 3D BSE for the qqq proton. Our goal is 
to write down ( and solve) the Master equation for three fermion quarks 
complete with all internal d.o.f.'s, in the presence of both qq and direct qqq 
forces which reads [S] : 



where the definitions for the various momenta, and the phase conventions for 
the quark propagators are those of [8], while the direct 3-quark interaction 
V qqq in the last term is new, and given by (1.2). The central problem is now 
the reduction of this Master equation (2.3) through the three steps (A, B, 
C) vide Sect (1.2), as outlined below. 



Piz]PiO = 



Pi = {Pi±,Piz} 



(2.2) 
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2.2 Reduction of Master Eq from ^ to $ 

Since permutation symmetry plays a crucial role for a qqq hadron for fermion 
quarks, we define at the outset a pair of internal variables (£ ; 77) with, say, 
the index #3 as basis, as [22] 

V / 2£ 3 =Pi-p 2 ; v 7 ^ = -2p 3 + Pi +p 2 ; ^ = ^1+^2+^3 (2.4) 

where the time-like and space- like parts of each are given by (2.2), and 
the corresponding 3-vector defined as pi = {pi±,Pi Z }. Two identical sets of 
momentum pairs £1,771 and £2,^2 are similarly defined, but can be expressed 
in terms of the set (2.4) via permutation symmetry. Now Step A, which is 
designed to dispose of the fermion d.o.f.'s for a qqq system, consists in defining 
an auxiliary scalar function $ related to the actual BS wave function \1/ by 

m 

y = n l23 s F l(- Pl Mp t p 2 p 3 )w(P) (2.5) 

Here the quantity W(P) is independent of the internal momenta but includes 
the spin-cum-flavour wave functions Xi °f t ne 3 quarks involved (see [?] : 



W(P) = lx'<P' + x"<P"]/V2 (2-6) 

where 0', 0" are the standard flavour functions of mixed symmetry [25] [not to 
be confused with the 3D wave function !], and x' >x" are the corresponding 
relativistic spin functions. The latter may be defined either in terms of the 
quark # indices as in Eqs (1.2) or (2.3), or sometimes more conveniently in 
a common Dirac matrix space as 



W >; \x" >= [ 2M [i^;i%/V3}C/V2} Pl (gi [[1;7 5 7>(P)] Q (2.7) 

where the first factor is the /?7-element of a 4 x 4 matrix in the joint spin 
space of the quark #s 1, 2 [26], and the second factor the a element of a 4 
x 1 spinor in the spin space of quark # 3; C is a charge conjugation matrix 
with the properties [27] 

-% = C-\C; 75 = C-^C; 

and 7^ is the component of 7^ orthogonal to P^. Finally, the representations 
of the flavour functions 0', 0" satisfy the following relations in the "3" basis 
[28] 

< 0"|l;r( 3 >|0" > = < 0'|l;_If(3)|0' > 
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2.3 Gordon Reduction on V qq z & V qqq 

At this stage we indicate the effect of Gordon reduction on the pairwise 
kernels Vfofji) and the 3-body kernel V qqq , following the original treatment of 
[29] . (also quoted in [8] ). The strategy lies in a close scrutiny of Eqs. (2.3-2.5) 
with a view to eliminate the Dirac matrices in favour of the Pauli matrices 
o^y. To that end, we express the Dirac propagators Sp(Pi) as 

iSpipi) = (m q - inf®.pi)/Aii A; = m 2 q +p 2 i 

Now employing the notation 

V® = {m q -i 1 ®.p i )i 1 f (2.8) 

the result of Gordon reduction is expressed by the formula [291 E] 

V W.V^ = -M 1 2 2 -(g+g0 2 +^(g+g0-(^i+^)x(g-g0-4 ) 4 2) (9-g , ),(9-g , ) fe 

(2.9) 

where we have used a mixed vector-tensor notation for the various 3- vectors 
on the light-front in the sense of (2.2) and a short-hand notation q for q 12 . It 
is now possible to identify the pairwise kernel arising from the first group of 
terms on the rhs of Eq.(2. 3) after expressing it in LF notation of Eq.(2.2) 
as follows: 

V qq3 = V(Li'i) = glW.vW x x 1 (2.10) 

where we have included a net color factor of (—4/3) arising from a folding 
of the factor F 12 for each qq pair in a 3* state, into A 3 /2 for the spectator 
quark # 3 in a 3 state; the corresponding gluon propagator has been taken 
on the light-front, in the notation of (2.2). Similarly the 3-body kernel V qqq 
in Eq.(2.3) may be identified as a sum of 3 terms arising from a folding of 
(1.2) with the 3 Dirac factors (m q — ij^.pi) in the second term of (2.3). 

V m = £(-4/3)^ (1) .(fc 2 - h) x V^.V^ + {2} + {3}]/ '{kfkfkl} (2.11) 

123 

using the notation V® of (2.8). The net effect of the (antisymmetric) com- 
bination of the color matrices A is represented by the factor (—4/3) sitting 
in front; and ki = Pi — p\. The 'bosonized' form of Eq. (2.3) is now 

*< p ^» = gAk/ift v,(4 ' e) * w ' ri - P3) 
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Figure 1: (a) Pictorial view of 2- & 3- body interactions; (b) 'Mercedez-Benz' 
diagram for qqq-force 

in terms of the structures (2.10) and (2.11) for the 2- and 3-body kernels 
respectively. And the Gordon reduction formula (2.9 ) may be repeatedly 
used to used to further simplify the structures of V^jEq (2.10)] and V qqq [ Eq. 
(2.12)], which are collected as follows. 

2.3.1 Further Reduction of V qq s &; V qqq 

The 2-body and 3-body interactions are pictorially represented by fig 1(a) 
where the vertices of the triangle ABC stand for the position vectors (r 1; r 2 , r 3 ) 
of the 3 quarks in configuration space, while the centroid O has the position 
vector R = ((ri + r 2 + i"3)/3. The corresponding momenta of the gluon 
lines connecting them for pairwise qq interactions are precisely the quanti- 
ties qij — q[j that appear as arguments in Eq. (2.10), while the momenta 
associated with the gluon lines connecting all the 3 quarks together, a la the 
F-shaped fig. (lb), are the quantities ki corresponding to the vectors (ri — R). 
In the notation of Eq.(2.4), the former momenta are precisely £j, while the 
latter are the dual variables fji. [This interpretation will prove useful in the 
subsequent analysis]. 

The pairwise terms V qqi were earlier considered in detail [8]. The 'con- 
vective' terms are of the form (pi +P2)-{p[ +P2) which simplifies to — M\ 2 m 
— 4M 2 /9, plus smaller terms. The spin-orbit terms are unimportant for the 
ground state of the qqq proton, but the spin-spin terms survive for its ground 
state S — 1/2; their angular average works out as : 

ij 
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where we have introduced the total spin operator E for the proton state, with 
S 2 = 3. Collecting these results the sum of the 3 pairwise interactions (2.10) 
work out as 

V qq s = 9^ x £ ' M2P " - 3p 32 (S 2 - 9)] (2.13) 

We now consider the more interesting term V qqq , Eq.(2.12). The numerator 
can be broken up into 3 distinct groups : a) convective; b) spin-orbit; c) spin- 
sp[in-spin. Using the notation 2pi = Pi+p'i, the convective part vanishes (due 
to antisymmetry with too few d.o.f.'s) : 

E4[pi.p 2 ][2p 3 -(A:i-A; 2 )]^0 

The spin-orbit part comes from two groups of terms, one of which vanishes: 

^22p 1 .(k 1 - k 2 )[a 3 .(p 3 x k 3 ) 

where a mixed notation is used for the Pauli matrices, light-faced tensor for 
the 4D and hatted vector for the 3D forms. The other group gives rise to a 
non-zero contribution as 

J2^Pi-P2[^hu(k^ - k 2tx )\ » (-4*M 2 /9)5>^3 x m (2.14) 

i 

in the notation of Eq.(2.4). Note the overall antisymmetry of the last factor 
(independent of the index 3). 

Finally the (totally antisymmetric) spin-spin-spin contribution SSS is 
defined as : 

SSS = ^[icr^kivihti - k 3 ^xpk 2p i(TxIha (2.15) 

123 

The simplification of this expression (which is elementary) is helped by the 
identity afy = e^a^ . The next step is angular averaging over each of the 
independent 3-vectors involved in each of the three terms above, namely, 
< kikj >= k 5ij/3. The resulting expression for SSS is 

sss = ^{kikiki}^ x^Ei 

Substituting for SSS in Eq (2.11) and ignoring the spin-orbit terms (for the 
ground state of the proton) the net 3-body force is 

K^^^x^E^ (2.16) 
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3 3D Reduction By Green's Function Method 

We are now in a position to implement Step C of our program ( establishing a 
3D-4D interconnection between the corresponding wave functions) , for which 
it is convenient to employ the Green's function approach [22]. Calling the 
4D Green's functions associated with ^ and $ by Gp and Gs respectively, 
the connection between them analogously to Eq.(2.5) may be written as 



G F (£w ?rf) = W(P) ® U 123 S^(- Pi )G s (^ ^ / )ni 2 3^ 1 (-^)W(P / ) (3.1) 



where we have indicated the 4-momentum arguments of the Green's functions 
involved, in a common 5*3 basis (£, rj), and expressed the spin-flavour depen- 
dence of Gf as a matrix product implied by the notation W(P) <8> W(P'). 

Now to show the 3D-4D interconnection it is enough to work at the level 
of the 'scalar' Green's function Gs ( relabelled as G for simplicity), since it is 
trivial to include the spin-flavour d.o.f.'s in a matrix notation via (2.5) later. 
The main steps for the scalar qqq Green's function are shown next, while 
Appendix A sketches the corresponding process for a typical qq subsystem, 
from 3D reduction to 4D reconstruction, so as to serve as a simpler prototype 
for the actual qqq case. 



3.1 Reduction from 4D G(£r);grf) to 3D G{£fj;£'fj') 



We shall now outline a Green's function method to establish a 3D-4D inter- 
linkage between the 4D $ and the 3D </>, following a generalization of the 
procedure developed some time ago [22], so as to include the effect of the 
direct 3-body force V qqq . To that end we shall rederive the qq denominator 
functions Dij of [22] arising from single integrations over d&o, to show that 
they are exactly proportional to the 3D denominator function -D123 arising 
from a double integration over d^ drj in the direct 3-body term to be consid- 
ered here. [This is a big improvement over the earlier derivation [22] where 
this property was missing, and helps pave the way for putting together the 
effects of both 2- and 3-body forces within a common dynamical framework 

]• 

The full 4D Green's function for 'scalar' quarks is G(^rj; £'rf) (taking out 
the cm 5-fn), while its 3D counterpart G is [22] 




(3.2) 
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where the time-like subscript '0' should be read as the corresponding LF 
component in the sense of Eq.(2.2). Note that both G and G are S3 sym- 
metric, since the combination d^odrjo has this property. But the two hybrid 
Green's functions defined as [22] 

Gsnfofk&iQ = J drtodrfnGitmM) (3.3) 

are not S3 symmetric (hence they are indexed), since the integration now 
involves only one of the two C,,i] variables. Next, the 4D G(£rj; £'rj') satisfies 
a BSE of the form 

i(2^G(^,ev') = E/z^-^g)^^;^) 

123 J 4A 1 A 2 

+ / (3-4) 

which is analogous to Eq. (2.12) for the corresponding wave functions, except 
for a change of variables wherein the factors 1/4 and 1/9 in the first and 
second groups of terms stem from the relations (2.4) among the corresponding 
variables. It is understood that both these types of interaction have covariant 
3D support, and they incorporate the effect of 'Gordon reduction' as outlined 
in Sect 2.3 above. The 3D reduction of (3.4) is now achieved by integrating 
it a la (3.2), leading to a structure of the form : 

(2vr) 3 G(^V) = E^r/^C&.aGCe^g^) 

123 ^ V ^12 J 

To explain the structure of certain factors, Eq (3.5) shows that we have two 
types of 3D denominator functions and -D123, associated with pairwise 
2-body and direct 3-body forces respectively. It is easily shown that they are 
simply related to each other. To that end, Appendix A already shows the 
structure of Dij : 

D12 = — p 12+ ; £>i2+ = 2M[ujf ± p 2+ + Ldj ± p 1+ - P12-P1+P2+] (3.6) 
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where cuf ± equals w? q + pj±- Using the on-shellness (A 3 = 0) of the spectator 
(#3) then gives P\ 2 - = P- ~Pz- which reduces further to P_ — uJ 3J _/p 3+ . Its 
substitution back in (3.6) leads to the S3 symmetric result: 

p 3+ D 12+ = T^D 123 = 2 Y,{P2+P3+^1±} ~ 2Pl+P2+P3+P- (3.7) 

Here we have anticipated the structure of -D123 associated with the V qqq term 
whose formal definition is 

Pldq 12 -dp, 



D 123 J 4M 2 (2i7r) 2 A 1 A 2 A 3 1 ' ' 

where the (double) contour integrations in the indicated variables leads to 
the desired result. The resultant 3D BSE for G is now 

1 OQ \/ J 



+ 



^-y 3 J drd^AtsVs^Ws) (3.9) 



The structure of this equation reveals some interesting symmetries, when 
(2.13) and (2.16) are substituted for V qq3 and V qqq respectively. Namely, in 
the first group of terms the i] variable is the spectator since the integration is 
only over the £ variable, while in the second group of terms, their relative roles 
are interchanged with the £ variable effectively a spectator ! This fact is not 
immediately apparent because of the double integration involved, but a little 
reflection shows that the absence of the £ variable in (2.16) effectively ensures 
its spectator status in the argument of the corresponding Green's function. 
We shall see this feature more clearly in the coordinate space representation 
of the equation (corresponding to Eq (3.9) for G, to be considered in Section 
5. But first we turn our attention to the reconstruction of the 4D Green's 
function from the 3D form G satisfying (3.9). To that end we note that 
the qq and qqq group of terms in (3.9) need different strategies and are best 
handled separately, one at a time, (temporarily) ignoring the presence of the 
other. 
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4 Reconstruction of 4D Wave Function 



4.1 Reconstruction of G with qq Forces Only 

Considering first the qq group, we proceed exactly as in (A. 10) of Appendix 
A for the qq system. Namely, first express G 3r] , eq.(3.3), in terms of the 3D 
G, using the notation of (2.2) and the result of (3.7): 

G 3v (^ = o- a i'v') . D> ™ A , (4.1) 
2«7rp 32 AiA 2 2iTrp' 3z A' 1 A 2 

In a similar way the fully 4D G function is expressible in terms of the hybrid 
function G 3 £ as 

<v> - g sfe;^ &«as&4 (4 ' 2) 

Now since the function is not determined from qqq dynamics alone, we 
invoke an ansatz similar to, but more symmetrical than, [22J : 

Gsc(&W iWs) = G(ifj; £'v')F(p 3 ,p 3 ) (4.3) 

where the balance of the p 3 (spectator) dependence is in the (as yet unknown) 
F function, subject to an explicit self-consistency check. To that end, try 
the ( symmetrical) Lorentz-invariant form 

F(p 3 ,p' 3 ) = ^5(A 3 -A 3 ) (4.4) 

and integrate both sides of (4.4) w.r.t. dp 30 dp 30 , in the notation of (2.2), 
to show that the consistency check is met with the mass shell value of the 
spectator momentum : 

A = fPfePL 

2m 

Substitution from A 3 in (4.4) then gives the symmetrical form 



F (P3, V'z) = 4p3 z p' 3 



, 5(A 3 - A> 3 ) 



2niA 



3 



whence the 4D G-fn in terms of G via the sequence (4.3) and (4.2), works 
out as 

123 A i A 2 A X A 2 {2mfJA 3 A' 3 
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4.2 Reconstruction of G with qqq Forces Only 

Next we consider only the qqq group of terms for a reconstruction from G ( 
Eq (3.5)) to G (Eq (3.4)). This case is more akin to the qq case of Appendix 
A in the sense that if we use the collective indices (£, n) = p and (£', n') = p' 
for the initial and final state arguments of the G-function we can define two 
kinds of G as follows. 

G(p; p>) = J d Po G(p; p'); G(p; p') = J dp' G(p; p') (4.6) 

where the integrals on the RHS are each of the double integral type, so that 
the definitions of the denominator functions involved are given by Eq.(3.8). 
Thus, following (A. 9), G and G are connected as follows 

together with a second one with the roles of p, p' interchanged. Continuing 
exactly as in Appendix A, G of Eq (3.4) gets expressed in terms of G on 
the RHS, since the interaction V qqq does not involve the variables p' . Thence 
another application of Eq.(4.7) for expressing G(p; p') in terms of G finally 
yields the desired connection: 

G(P ' P } " (2i7r) 2 A 1 A 2 A 3 P A[A' 2 A' 3 (48) 

where the symbol p stands collectively for (£, rj), as does p'. Note that in this 
case the reconstruction of G in terms of of (the fully reduced) G is unique, 
and does not require any extra ansatz like (4.4). 

4.3 Putting Both V qq & V qqq Together 

We have now two distinct types of 3D-4D interconnections valid for pure qq 
and pure qqq forces respectively. But when both types of forces are present, 
an exact interconnection is very difficult to derive. We shall therefore strive 
for an approximate but sufficiently realistic solution based on the relative 
strengths of the two forces, namely, an overwhelming preponderance of qq 
over qqq forces. To give effect to this strategy, we first note that the 3D 
level does not involve any approximation since Eq. (3.9) treats both types 
of interaction on par. It is only at the 4D level of reconstruction that an 
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approximation is necessary for putting the two forces together. Now taking 
account of the dominance of qq over qqq forces, the simplest choice is to 
prefer the structure of Eq (4.5) over that of Eq. (4.8). And although the 
ratio of the two 3D-4D interconversion jackets in (4.5) over (4.8) is a singular 
quantity 

6(A 3 -A' 3 ) x (2ttz) v /a^ 

it causes no harm for physical amplitudes since such singular quantities get 
ironed out on integration over the various momenta [23] . With this un- 
derstanding, Eq.(4.5) represents the reconstruction of the full 4D Green's 
function G(£rj; £'//) m terms of the 3D quantity G(^fj; £'77') where the latter 
now satisfies Eq (3.9) which includes the effect of both qq and qqq forces. 



4.4 3D-4D Interconnection for Wave Functions 

Finally the spectral representations of G(£r]; £'t/) an d G(£f); for the qqq 
system, exactly on the lines of (A. 12- A. 13) for a qq subsystem near a bound 
state pole P 2 = —M 2 , are 

G{tm i'rf) = E i)/( p2 + M2 ); ( 4 - 9 ) 

n 

and a similar equation for G vis-a-vis 0. These give the connection between 
the 4D wave function <3> which satisfies Eq (2.12), and the 3D wave function 
corresponding to the Green's function G which satisfies (3.9). For purposes 
of evaluating transition amplitudes via Feynman diagrams, it is convenient 
to index $ as $i + $2 + $3, and a corresponding indexing for the associated 
vertices as V = V\ + V 2 + V 3 , as in ref [22], so as to keep track of which 
quark is involved in which vertex. 



rh ft \- ^3 gl23 .ft ^ yjS(A 3 ) 

AiA 2 A 3 A X A 2 (2m) b ^y/ A 3 

The 5-function in eq.(4.7) has nothing to do with any connectedness problem; 
see ref. [22] for detailed reasons. Finally, the use of Eq.(2.5) with (4.10) yields 
an explicit structure for the actual (fermionic) wave function \I/ as 



123 \£K%) 
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5 Complete cj) Equation In Coordinate Space 

Our final task is to set up (and solve) the 3D BSE for which may be inferred 
from (3.9) by making use of a spectral representation similar to (4.9), and 
going to the pole P 2 = -M 2 : 

(27r) 3 ZW(^) = E^= / dH'X^dlvs) 

123 V ^ J 

+ (5 ' 1) 

where V qq3 is given by (2.13); V qqq by (2.16); and D 123 by (3.7). To transform 
this equation in coordinate space, define the combinations analogous to (2.4) 

as 

v/2s 3 = n - r 2 ; V6t 3 = -2r 3 + r l +r 2 (5.2) 

Then a Fourier transform to the coordinate — space representation gives for 
the pairwise terms V qq3 the coulombic structure 

y » (4) = / H exp [tg.s 3 ] x E Jgdg. _ | (S . _ 9)p3 J (5 ,) 

which multiplies the coordinate space wave function 0(s, t) ( S3 symmetric 
in its arguments). Similarly, the (S3 symmetric) double Fourier transform of 
V qqq is defined as 

v ' 3 ^ ( > = I Sgr exp K '' s + itlA x § CT [ *" ( * 2 x ^ 

(5.4) 

The integration in (5.3) involves a 3D 5-function <5 3 (s3) ; and the double 
integration in (5.4) additionally involves the factot l/l^l due to the 77 in- 
tegration. The 5 function being singular may be 'regularized' by using a 
differential representation (acting on the scalar function 0(s, t)), with \s 3 \ = 
N : 

The result of integration in (5.3) is then expressed as 



£<\S 3 \ ^' P3 



(5.6) 
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Next the 3-body term term (5.4) integrates, with the help of (5.5), to 



Now the product 1 / s^t^ can be approximately replaced by the S3 symmetric 
expression 2/[s| + t 2 } which can be taken out of the summation sign (with 
index '3' dropped) to give ^2d 2 3 which in turn equals (3/2)(<9 2 + <9 2 ), an S3 
symmetric sum with index '3' dropped again. This finally leads from (5.4) 
to the S3 symmetric form of the 3-body term in coordinate space: 

For further manipulations it is convenient to replace the antisymmetric spin 
operator A = [i<7i.(<7 2 x 63)] by one of its possible eigenvalues as follows: 
Squaring this quantity yields in a simple way 

A 2 = -A - 15 + £ 2 ; S = o-i + o- 2 + 6- 3 

from which one obtains the successive results 

A 3 + A 2 + 15A = £ 2 A; A 2 (A + l) 2 = [£ 2 - 15] 2 

The last equation yields the four solutions 

A =~\± y/l5-&-l/4; A = ~ ± v'S 2 - 15 - 1/4 (5.8) 

whose substitution in Eqs (5.4) and (5.7) summarises the full content of the 
total spin effect of the V qqq term. Finally, the denominator term D123 is a 
differential operator in coordinate space : 

£123/4 = £[(-™? + dl ± )d lz d 2z ] - tMd lz d 2z d 3z (5.9) 

123 

Thus satisfies the following equation in coordinate space : 

£i230M) = [V (2 \s) + V^ 3 \s,t)}<p(s,t) (5.10) 
where the various operators are given by (5.6) - (5.9). 
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5.1 Simplified Form of the Equation (5.10) 

Since this paper is intended as a preliminary mathematical framework for 
the dynamical effect of the spin terms, especially the spin-rich V qqq term 
on the spatial structure of the proton wave function, pending a full-fledged 
study of the proton spin anomaly, we shall at this stage merely outline a 
qualitative procedure to determine the nature of the solution of Eq. (5.10), 
with particular reference to the possible role of the eigenvalues (5.8) of A, on 
the solution of this differential equation. To that end, we shall make some 
drastic simplifications, starting with the differential operator (5.9) whose 
principal terms (up to second order ), in momentum space, work out as 

9x r 2M 2 9l/ . 9N AM 2 , 9 
^123 « — (£ + Vl) + [~^-~ 2m 2 Ml + r>l) + - M 79) 
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Further simplification arises with the 'special' value M = 3m q which then 
yields a simple yet transparent expression in coordinate space : 

AM 2 

Dus * —l-d 2 s - d?] (5.ii) 

an operator with full rotational invariance and S3 symmetry in the 3D space 
( on light-front) defined by Eq. (2.2). In a similar vein, the 2-body terms 
(5.6) can be simplified by the replacements p iz ~ M/3, so as to yield 

W(s, t) = 8 ^ M3a ° + ^(f-^ [fl g + dn (5.12) 

where we have made a further simplification based on certain standard in- 
equalities [31] 

And the 3-body term (5.7) may be written more compactly in terms of of 
the operator A with eigenvalues (5.8) as 

One has now a differential equation for (5.10), with the simplified operators 
(5.11-13) which exhibit a 6D symmetry. Taking R 2 = s 2 + 1 2 , the 6D Lapla- 
cian for the ground state of the proton (with all angular d.o.f.'s dropped) 
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takes the simpler form 



[4M /9 + 27*/2 + 5^ l<9 « + R dR)m) + = ° 

which on rescaling with the dimensionless variable X = MR leads to the 
simpler form 

5.2 Spin Effects of V qq & Terms on qb Singularity 

The first thing to notice from the Equation, (5.14), is that the spin- 
dependent parts of both V qq and V qqq appear in the multiplying factor with 
the 6D Laplacian acting on 0, and are proportional to a s and a 2 respectively, 
with the V qq term having a milder singularity (~ X~ r ) than the V qqq term 
(~ X~ 2 ). Further information on the singularity of the differential equa- 
tion (at points other than R = 0) hinges on the nature of the eigenvalues of 
the spin operator A which appears in above multiplying factor. Now these 
eigenvalues are given by (5.8), two of which are complex for the state of the 
proton (£ 2 = 3), but the other two are real. Of the real solutions the positive 
eigenvalue does not yield a zero in this multiplying factor but the negative 
eigenvalue (~ —4) does give a zero for a real value of R at a point Xq between 
X 2 = and X 2 = oo : 

X 2 - -^X + = 0;X = ± ] J[a 2 /32 - (5.15) 

where the value of S 2 = 3 for the proton state has been substituted. Indeed 
this zero ( for A ps —4) is a key element of the dynamical effect of the spin 
- rich 3-body force term, although the spin effect of the 2-body V qq term is 
marginal. 

To study this effect more quantitatively, we seek an approximate solution 
of Eq. (5.14) in the neighbourhood of the zero at (5.15) for which a crude nu- 
merical estimate suggests the following location. Taking a 3-flavour structure 
a s = 27r/[91nM/A], with A ss 150MeV, one finds a s w 0.39, whence 

X w 0.0689 ± Xi; X 1 = 0.5841 
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This shows that the qq term gives ~ 10% shift around the central value of 
Xi, which corresponds to R ~ 0.12/m, and is almost entirely due to the spin 
effect of the V qqq term. Therefore in the spirit of this qualitative investigation, 
it makes sense to drop this 10% effect, in which case Eq (5.14) simplifies to 

+ ^ + ^-0; X 1 = 0.584 (5.16) 

To bring it nearer to a standard form, transform to the independent variable 
z according to zXq = Xq — X 2 , which yields 

a 



z(l - z)dr z (f) - 3zd z (f> - /3VT^z4> = 0; (3 = -^=X 1 « 0.058 (5.17) 

3v 2 

This equation is almost (not quite) of the hypergeometric form but it can be 
reduced to a standard one (with singularities located at z = 0, 1, oo) [32] by 
exploiting the smallness ( (3 = 0.058 ) of the last term to replace it with a 
constant ( with value corresponding to z = 1/2) : 

«(1 - z)<9^ - 3zd z <p - 4=0 = (5.18) 

v2 

An equivalent equation may be obtained with the transformation z = 1 — x : 

ar(l - x)fi£0 + 3(1 - x)d x <j> - -^=(p = (5.19) 

V2 

where x = 1 corresponds to the location X = X\ of the singularity, as 
indicated in Eq. (5.16). [ Note that this singularity corresponds to the 
negative eigenvalue of A, thus reflecting the dynamical effect of the spin-rich 
V qqq term, while a positive or complex eigenvalue of A would result in the 
disappearance of this singularity ]. The solution of Eq.(5.19) is then given 
by [32] 

= F(o,6|3|ar); a + b = 2; ab = (5.20) 

v2 

with the entire machinery of hypergeometric functions available [32] for ex- 
ploring its properties according to need. The scaled variable x is related to 
the 6D distance R by MR — X\x where X% w 1.5a s is directly related to 
the location of the singularity induced by the spin structure of the Y-shaped 
3-body force . And the reconstruction of the full Bethe Salpeter wave func- 
tion \I/ as a sum of 3 pieces \&j is now only a matter of substituting (5.20) in 
(4.11), whence the corresponding vertex functions Vi are immediately iden- 
tified [221 [23] for specific transition amplitudes a la Feynman diagrams. 
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6 Retrospect : Critique Of The 3-Body Force 



In retrospect, we have considered, in conjunction with pairwise qq forces, the 
effect of a new type of 3-body force V qqq , rich in spin content, on the ana- 
lytical structure of the qqq wave function in the high momentum regime of 
QCD where the confining interaction is unimportant, rendering the dominant 
force Coulombic. As to the anatomy of this (spin-rich) V qqq , we have taken 
it to be generated by a ggg vertex ( a genuine part of the QCD Lagrangian 
) wherein the 3 radiating gluon lines end on as many quark lines, giving rise 
to a (Mercedes-Benz type) ^-shaped diagram, a la fig 1. From a physical 
point of view it is natural to expect that such a spin-rich structure should 
play a potentially crucial role in the so-called 'spin anomaly' of the proton, 
a subject that seems once again to have raised its head in the context of new 
polarized beam techniques now available [5] for resolving the issue experi- 
mentally. With that end in view, our strategy has been to determine the 
dynamical effect of the spin dependence of V qq and V qqq forces on the ana- 
lytical structure of the internal 3D wave function 0. [It is emphasized that 
this effect must be carefully distinguished from the " kinematical " effect of 
spin, which manifests in several other ways, namely through the presence of 
various 7 matrices that appear in equations like (4.11) connecting <p to the 
full 4D wave function *&]. Indeed we found in Section 5 that while the spin- 
dynamical effect of 2-body forces is marginal, that of the spin-rich 3-body 
force is quite pronounced, and shows up through the possible " eigenvalues 
", Eq (5.7), of the spin operator {io\.Oi x cr 3 ) which is a part of V qqq : only 
a negative eigenvalue (there is only one !), induces a singularity in the dif- 
ferential equation for 0, but not others. The resulting dynamical effect is 
expressed by a hypergeometric function, Eq.(5.20), which gets folded into 
the 4D BS amplitude ^/ via Eq.(4.11), thus fulfilling the original (limited) 
motivation behind this study. 

Unfortunately, the dynamical framework underlying the contents of this 
paper has had a long history born out of the author's long involvement with 
the so-called Bethe-Salpeter Equation (BSE), often with extended periods of 
stalemate (and frustration !), but mostly centred around the quest for a sat- 
isfactory definition of probability within the BSE framework. Eventually it 
became possible to settle for a toned down version of BSE, that of a Salpeter- 
like equation ( 3D support for the kernel), which is amenable to a probability 
interpretation at the 3D level. And its 4D features, although present in the 
original formulation itself [6], had for decades remained hidden from view, but 
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were finally dug out [7] in the context of an independent approach designed to 
explore a 3D-4D interconnection between the corresponding BS amplitudes 
when the kernel has a 3D support [8]. Further, the lack of covariance in the 
original formulation [6J was subsequently remedied via a special (instanta- 
neous) frame of reference in which the composite hadron of 4-momentum P M 
is at rest [33]. This result turned out to be in conformity with the Markov- 
Yukawa Transversality Principle (MYTP) [TTJ [T2], which happens to be a 
'gauge principle' in disguise [13] : It ensures that the interactions among the 
constituents be transverse to the direction of P M . Subsequent refinements 
have been mostly technical , especially the use of Dirac's light-front formula- 
tion so as to extend the dynamical range of validity of the BS framework by 
overcoming certain practical problems like 'Lorentz mismatch disease' [T7] 
associated with different vertices of a given Feynman diagram. 

It therefore looked worthwhile to employ this old-fashioned formalism 
in the present context of a new kind of 3-body force V qqq , with enough de- 
tails put in for a reasonably self-contained description without having to dig 
frequently into the original sources. This has also given an opportunity to 
make several refinements, especially the derivation of a common denomina- 
tor function associated with both the qq and qqq forces and extending the 
earlier formalism to accommodate the V qqq force, for which this theoretical 
framework seems to be well suited. 

The next part of this programme involves the derivation of the requi- 
site QCD parameters for the baryon spin, for which some key ingredients 
are i) the forward scattering amplitude off the £7^75 operator, inserted at 
the individual quark lines, and ii) a more elaborate 2-gluon contribution to 
the proton spin, which is reserved for a subsequent paper within the same 
formalism. 
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Appendix A : qq Subsystem Formalism 

Using the correspondence (2.2) of the text, most of the covariant instant form 
results of [22] may be taken over to the present light-front situation |7j. In this 
Appendix, we outline the structure of the 4D / 3D interconnection between 
the corresponding Green's functions for the qq sub-system as a prototype for 
the actual qqq system considered in Section 2 of the text. The qq Green's 
functions satisfy the respective equations [22J: 

(2n)HG(q, q'; P) = J d 4 q"V(q, q")G(q", q'; P); (A.l) 

G(q,q')= Jdq dq' G(q,q';P) (A.2) 

Integrating both sides of (A.l) gives via (A.2), the 3D BSE for a bound state 
(no inhomogeneous term !): 

(2n) 3 D(q)G(q, q') = J d*q"V(q, q")G(q", q') (A.3) 

where the 3D denominator function D(q) is defined as 

2m r dqo (A.4) 



D(q) J AxA 2 
leading (for general unequal mass kinematics) to [8] 

*«> = D + (V = 2P4f- KM )Zl' m ^ (A- 5 ) 

where A is the triangle function of its arguments. Now define the hybrid 
Green's functions 1221: 



G(q, q') = J dq G(q, q'; P); G(q, q') = J dq' G( qi q>; P) (A.6) 
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Using (A. 6) on the RHS of (A.l) gives 

(2n)HG(q, q'; P) = — L- J d 3 q"V(q, q")G(q", q') (A.7) 
Integrating (A.l) w.r.t. dq' only, and using (2.7) again, gives 

(2n)HG(q, q') = — L- [ d 3 q"V(q, q")G(q", q') (A.8) 



A1A2 

From (A.8) and (A. 3), G and G are connected as: 



Interchanging q and q' in (A. 9) gives the dual result 

d ^=2^k diM) < A ' 10) 

Substitution in (A.7) gives the desired 3D-4D interconnection 

Next, spectral representations for the 4D and 3D G-fns in (A. 11) give [22] 
G(q, q'; P) = E ®n(q; P)KW\ P )/(P 2 + M 2 ); (A.12) 

n 

G(q,q') =J2Mq)€m/(P 2 + M 2 ) (A.13) 

n 

where and <p n are 4D and 3D wave functions, so that their interconnection 
(valid for any bound state n), is expressed by: 

r W ).A 1 A 2 $(,;P) = ^|M (A.14) 

which tells us that the covariant vertex function T on the light-front is a 
function of q only [22] via its definition (2.2). This derivation is a prototype 
for the qqq case of text. 
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